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The  reaction  between  Hj  and  O2  has  been  studied  in  a  reflected  shock  tube  apparatus  between  tem¬ 
peratures  of  1662-2097  K  and  pressures  of  400-570  torr  with  Kr  as  the  diluent  gas.  O  atom  atomic 
resonance  absorption  spectrometry  (ARAS)  was  used  to  observe  absolute  [O],  under  conditions  of  low 
[Hjlo  so  that  most  secondary  reactions  were  negligible.  Hence,  the  observed  [O],  was  the  direct  result  of 
the  rate  controlling  reaction  between  Hg  and  O2.  Three  different  reactions  were  considered,  but  experi¬ 
mental  and  ah  initio  theoretical  results  both  indicated  that  the  process,  H2  +  O2  — >  H  -t  HO2,  is  the 
most  probable  reaction.  After  rapid  HO2  dissociation,  O  atoms  are  then  instantaneously  produced  by  H 
-t  O2  ^  O  -t  OH.  Using  the  ab  initio  result,  conventional  transition  state  theoretical  calculations  (CTST) 
with  tunneling  corrections  give  the  expression  k*  =  1.228  X  10“^®  exp( —  26,926  K/T)  cm^  mol- 
ecule^^  s“^,  applicable  between  400  and  2300  K.  This  theoretical  result  agrees  with  the  present  experi¬ 
mental  determinations  and  those  at  lower  temperature,  derived  from  earlier  work  on  the  reverse  reaction. 


Introduction 

The  branching  chain  oxidation  of  H2  with  O2  is 
one  of  the  most  studied  reactions  in  combustion.  It 
was  recognized  very  early  that  branching  could  only 
occur  after  some  initiation  process  produces  radical 
chain  centers  [1].  Over  the  years,  five  initiation  re¬ 
actions  have  been  considered: 


H2  +  02 

^  H  -k  HO2 

(1) 

Ha  +  O2  ^  OH  -k  OH 

(2) 

H2  +  O2  “ 

^  0  +  H2O 

(3) 

Ha  +  M^H-kH-kM 

(4) 

O2  H“M~*’0  +  0  +  M 

(5) 

Reactions  1  and  2  have  been  the  preferred  initiation 
processes  since  reaction  3  requires  multiple  bond 
rearrangement  and  both  diatomic  dissociation  reac¬ 
tions;  4  and  5  require  substantially  higher  energy  [2] . 
Reaction  2  was  almost  always  the  choice  in  shock 
tube  induction  delay  experiments  [3-6]  until  about 
1985,  with  the  value  from  Ripley  and  Gardiner  [3], 
k2  =  4.16  X  10“^^  exp(  — 19,605  K/T)  cm^  mole¬ 
cule”^  s”^,  being  typical.  However,  later  flow  tube 
studies  at  room  temperature  [7,8]  on  the  reverse  of 
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reaction  1,  H  -k  HO2,  ^  H2  +  O2,  gave  a  rate 
constant  of  A:  _  1  =  6.7  X  10”^^  cm^  molecule”^  s”^ 
with  at  least  an  error  of  ±  50%,  essentially  corrob¬ 
orating  earlier  values  [9,10].  Using  Janaf  calculated 
equilibrium  constants  [11],  Warnatz  [12]  then  esti¬ 
mated  rate  constants  for  reaction  1  and  suggested  its 
importance  as  an  initiation  process.  With  derived 
773  K  results  from  Baldwin  et  al.  [13]  (i.e.,  A_i  = 
4.65  X  10”^^  cm^  molecule”^  s”^),  Tsang  and 
Hampson  [14]  further  refined  the  estimate,  albeit 
with  a  large  uncertainty,  suggesting  ki  =  2.4  X 
10”^°  exp(  — 28,500  K/T)  cm®  molecule”^  s”®  for 
300-800  K. 

From  the  above  discussion,  it  is  clear  that  cur¬ 
rently  used  values  for  initiation  with  reactions  1,  2, 
or  3  are  based  on  a  relatively  small  experimental  da¬ 
tabase.  In  the  present  work,  we  studied  the  H2  + 
O2  reaction  with  atomic  resonance  absorption  spec¬ 
trometry  (ARAS)  under  low  [H2]  conditions  using 
the  reflected  shock  tube  technique.  Absolute  rates 
of  O  atom  formation  are  measured  during  the  initial 
stages  of  reaction  under  conditions  where  the  rate 
of  initiation  is  rate  controlling  and  nearly  chemically 
isolated.  Previous  ARAS  studies  [6,15,16]  did  not 
use  a  sufficiently  high  [02]/[H2]  ratio  to  allow  kinetic 
isolation  of  the  initiation  process.  The  difficulty  in 
assessing  initiation  with  ignitable  quantities  of  H2 
and  O2  is  well  illustrated  in  the  work  of  Pamidimuk- 
kala  and  Skinner  [6],  who  finally  concluded  that  ther¬ 
mal  dissociation  of  impurities  probably  initiated  the 
branching  chain  oxidation  in  their  work.  By  contrast, 
in  the  present  work  [H2]  is  so  small  and  impurity 
levels  are  so  low  that  any  observed  atom  formation 
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ELEMENTARY  REACTION  KINETICS 


must  derive  from  a  direct  bimolecular  reaction  be¬ 
tween  H2  and  O2. 


Experimental 

The  present  experiments  were  performed  with 
previously  described  equipment  [17],  and,  there¬ 
fore,  only  a  brief  description  of  the  system,  along 
with  those  features  unique  to  the  current  experi¬ 
mental  procedures,  will  be  presented  here. 

The  apparatus  consists  of  a  7  m  (4  in.  o.d.)  304 
stainless  steel  tube  separated  from  the  He  driver 
chamber  by  a  0.004  in.  unscored  1100-H18  alumi¬ 
num  diaphragm.  The  tube  was  routinely  pumped 
between  experiments  to  less  than  10“®  torr  (760  torr 
=  1.01325  X  10®  Pa)  by  an  Edwards  Vacuum  Prod¬ 
ucts  model  CRIOOP  packaged  pumping  system.  The 
velocity  of  the  shock  wave  was  measured  with  eight 
equally  spaced  pressure  transducers  (PCB  Piezo- 
tronics,  model  113A21)  mounted  along  the  down¬ 
stream  part  of  the  test  section  of  the  shock  tube  and 
recorded  with  a  4094C  Nicolet  digital  oscilloscope. 
Temperature  and  density  in  the  reflected  shock  wave 
regime  were  calculated  from  this  velocity.  This  pro¬ 
cedure,  which  includes  corrections  for  boundary 
layer  perturbations,  has  been  fully  described  else¬ 
where  [17,18].  The  digital  oscilloscope  was  triggered 
by  pulses  derived  from  the  last  velocity  gauge  signal. 
The  photometer  system  was  radially  located  6  cm 
from  the  end  plate.  All  optics  were  made  from 
MgF2.  The  resonance  lamp  beam  intensity  was  mea¬ 
sured  by  an  EMR  G14  solar  blind  photomultiplier 
tube  and  recorded  with  the  oscilloscope. 

The  technique  used  for  the  detection  of  the  tran¬ 
sient  O  atoms  is  atomic  resonance  absoiqition  spec¬ 
troscopy  (ARAS).  The  level  of  light  intensity  in  the 
source  is  so  low  that  photoinitiation  is  not  possible 
in  the  present  experiments.  In  earlier  work  [19],  an 
O  atom  curve  of  growth  was  determined  using 
Xo,  =  1  X  10  ®  in  1.8  torr  of  purified  grade  He 
at  50  W  microwave  power  to  give  an  effective  lamp 
temperature  of  490  K  [20] .  This  curve  of  growth  was 
slightly  modified  with  the  inclusion  of  additional  data 
obtained  in  a  later  study  [21].  The  present  experi¬ 
ments  were  carried  out  in  exactly  the  same  way  as 
in  this  earlier  study  [21]  except  that  H2  was  used  as 
the  reactant  instead  of  CH3.  Following  procedures 
discussed  earlier  [21],  28  kinetics  experiments  were 
carried  out  between  1662  and  2097  K. 

Gases 

High  purity  He  (99.995%),  used  as  the  driver  gas, 
was  from  Air  Products  and  Chemicals.  Scientific 
grade  Kr  (99.999%),  the  diluent  gas  in  reactant 
mixtures,  was  from  Spectra  Gases.  In  Kr,  the  ~10 
ppm  impurities  (N2,  2  ppm;  O2,  0.5  ppm;  Ar,  2  ppm; 
CO2,  0.5  ppm;  H2,  0.5  ppm;  CH4,  0.5  ppm;  H2O, 


0.5  ppm;  Xe,  5  ppm;  CF4,  0.5  ppm)  are  all  either 
inert  or  in  sufficiently  low  concentration  so  as  to  not 
perturb  O  atom  profiles.  Ultra  high  purity  grade  He 
(99.999%  from  AGA  Gases)  was  used  for  the  reso¬ 
nance  lamp.  High  purity  O2  (99.995%)  for  the 
atomic  filter  was  from  AGA  Gases.  Scientific  grade 
O2  (99.999%)  and  H2  (99.9999%)  for  reaction 
mixtures  were  obtained  from  MG  Industries  and 
were  used  without  additional  purification.  Test  gas 
mixtures  were  accurately  prepared  from  pressure 
measurements  using  a  Baratron  capacitance  mano¬ 
meter  and  were  stored  in  an  all-glass  vacuum  line. 


Results 

With  the  relatively  high  levels  of  O2  used  in  the 
present  experiments,  there  is  an  observable  absorp¬ 
tion  of  the  resonance  light  at  130  nm  by  molecular 
oxygen.  Since  this  absorption  is  uniform  over  the  O2 
bandwidth  and  [O2]  does  not  significantly  change 
during  an  experiment,  any  observed  decrease  in  the 
transmitted  light  intensity  without  reactant  H2  re¬ 
flects  an  increase  in  total  density  (due  to  vibrational 
relaxation  of  O2)  and/or  O  atom  formation  from  O2 
I-  M  ^  20  I-  M  at  high  temperature.  Hence,  we 
carried  out  a  limited  set  of  O2  dissociation  experi¬ 
ments  with  the  same  [O2]  that  was  used  in  the  ki¬ 
netics  experiments.  Under  conditions  where  O  atom 
formation  was  negligible,  density  relaxation  (gener¬ 
ally  <250  ,us  at  T  >  1600  K)  was  easily  determined. 
The  apparent  absorbance  due  to  the  vibrational  re¬ 
laxation  of  O2  was  a  small  correction  and  was  not 
strongly  dependent  on  temperature.  This  correction 
was  then  point-by-point  subtracted  from  the  higher 
temperature  runs  where  O  atoms  were  definitely 
formed  from  dissociation  as  reflected  in  a  substantial 
strong  increase  in  absorption.  The  corrected  absorb¬ 
ance  was  then  converted  to  [OJt  with  the  previously 
determined  curve  of  growth  [21].  Using  the  expres¬ 
sion  kj  =  Ro/2  [02][M],  where  Rq  is  the  rate  of 
formation  of  O  atoms  (i.e.,  the  slope  from  the  ex¬ 
perimental  [0](  against  t  plots),  values  for  were 
determined.  For  T  >  2200  K  the  results  were  within 
±  40%  of  those  accurately  determined  previously  by 
Jerig,  Thielen,  and  Roth  [22].  Hence,  we  have 
adopted  their  value  in  the  kinetics  model  for  the  H2/ 
O2  experiments.  In  the  kinetics  experiments,  we 
similarly  corrected  for  O2  density  relaxation  by 
point-by-point  subtraction.  In  general,  the  contri¬ 
bution  of  this  correction  was  ~20%;  that  is,  80%  of 
the  signal  was  due  to  chemical  reaction  with  reactant 
H2.  After  conversion  to  [0]t  using  the  curve  of 
growth  [21],  O  atom  concentration  profiles  were  de¬ 
termined,  two  of  which  are  shown  in  Fig.  1.  These 
experimental  profiles  were  then  compared  to  those 
generated  by  numerical  integration  of  the  reaction 
mechanism  listed  in  Table  1. 


RATE  CONSTANTS  FOR  Hj  +  O2  ^  H  +  HO; 
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Fig.  1.  Comparison  of  measured  and  simulated  O  atom 
profiles  for  initiation  by  reaction  1  (solid  line)  or  reaction 
2  (dashed  line).  Top  panel:  Pj  =  10.94  torr,  =  2.915, 
P5  =  2.549  X  10“  Tg  =  1981  K,  [H2]  =  6.357  X  10“, 
and  [O2]  =  1.561  X  10^^,  all  in  molecules  cm  =  1.6 
X  10“^®  cm^  molecule^^  s“^  or  k.2  =  7.3  X  10“^®  cm® 
molecule^®  s“®.  Bottom  panel:  Pi  =  10.94  torr,  = 
2.775,  ps  =  2.442  X  10®®,  Tg  =  1810  K,  [Ha]  =  1.227  X 
lO®®,  and  [O2]  =  1.508  X  10®®^,  all  in  molecules  cm“®,  ki 
=  8.0  X  10 “®^  cm®  molecule “®s^®  or  ita  =  2.25  X  10“®® 
cm®  molecule  “®  s“®. 


Discussion 

Experiment 

[OJt  profiles  (0  to  1  ms)  for  28  experiments  were 
generally  simulated  over  the  concentration  range,  0 
to  ~1  X  10®®  atoms  cm“®,  using  the  mechanism  of 
Table  1.  For  all  conditions,  the  only  important  re¬ 
actions  in  the  table  are  reactions  1,  2,  or  3,  possibly 
followed  by  reactions  6  through  9.  In  the  initial 
stages  of  reaction,  all  radical-radical  processes  are 
negligibly  slow  as  are  the  reactions,  H  and/or  O  -f 
H2O,  both  of  which  have  high  activation  energies. 
H2  dissociation,  reaction  4,  is  also  too  slow  for  the 
small  [H2]  employed.  However,  O2  dissociation  [22] 
can  significantly  contribute  to  O  atom  formation  at 
T  >  2100  K,  and,  therefore,  no  experiments  above 
this  temperature  were  included  in  the  database.  Un¬ 
der  all  experimental  conditions,  reactions  6  and  7  are 
effectively  instantaneous,  leaving  only  reactions  8 
and  9  contributing  to  the  time  dependence  of  [Olf. 


Hence,  fitting  is  dependent  only  on  the  values  of  the 
rate  constants  for  reactions  1, 2,  or  3  and  by  reactions 
8  and  9,  both  of  which  are  well  known  over  the  pres¬ 
ent  T  range  [25]. 

Even  though  most  of  the  reactions  are  relatively 
unimportant,  the  entire  set  of  reactions  in  Table  1 
was  used  to  simulate  [0]t.  Fig.  1  shows  two  such 
profiles.  The  solid  or  dashed  lines  are  simulations 
under  the  separate  assumptions  that  reactions  1,  2, 
or  3  initiate  the  reaction.  The  yield  of  O  atoms  per 
H  atom  destroyed  using  1  is  initially  2  since  both 
reactions  6  and  7  are  instantaneous.  Hence,  simu¬ 
lations  using  reaction  3,  where  only  one  O  atom  is 
formed,  are  identical  to  reaction  1  except  that  the 
rate  constant  for  reaction  3  has  to  be  two  times  that 
for  initiation  by  reaction  1.  In  all  cases,  the  values 
resulted  from  iterative  fits,  and  doubling  or  halving 
the  rate  constants  for  the  processes  results  in  values 
of  [0](  at  I  ms  that  are  about  twice  or  one-half  of 
those  from  the  best  fits.  Given  the  substantial  ex¬ 
perimental  error  in  absolute  [O]  (i.e.,  ~  ±  I5%)[21], 
the  fits  give  rate  constant  values  with  an  estimated 
accuracy  of  ~  ±  20%.  Inspection  of  Fig.  I  shows  that 
reactions  I,  2,  or  3  (or  a  combination)  could  account 
for  the  results,  and,  therefore,  we  determined  rate 
constant  values  for  the  entire  set  of  experiments  us¬ 
ing  each  of  the  three  reactions  as  the  initiation  pro¬ 
cess.  Starting  with  reaction  2,  if  this  reaction  is  pre¬ 
sumed  to  be  the  initiation  process  then  the  derived 
rate  constants  follow  the  Arrhenius  expression  ^2  = 
l.OI  X  I0“®  exp(  — 33,056  K/T)  cm®  molecule”® 
s”®,  whereas,  with  reaction  I,  the  value  is 

ki  =  1.50  X  10"®®  exp( -22,905  K/T) 

cm®  molecule"®  s"®  (6) 

for  1662  <  T  £  2097  K.  Considering  reaction  3,  if 
this  were  the  initiation  process  then  the  only  impor¬ 
tant  reaction  would  be  reaction  3  alone,  and,  in  this 
case,  the  values  for  would  be  ~2  times  equation 
6.  Since  the  data  are  obtained  over  such  a  limited  T 
range,  deciding  between  reactions  1,  2,  and  3  is  im¬ 
possible  from  just  these  results  alone.  A  choice  might 
be  possible  by  examining  A  factors  and  activation 
energies  from  the  three  analyses.  If  reactions  2  or  3 
were  operative,  then  one  would  expect  low  A  factors 
(i.e.,  ~I0"®^  cm®  molecule"®  s"  )  for  four-center 
or  tight  transition  states,  respectively  [27],  but  much 
higher  values,  A2  ~  I  X  10"®  or  A3  ~  3  X  10"®® 
cm®  molecule"®  s"®,  are  required  to  explain  the 
data.  Moreover,  it  is  also  worth  noting  that  reaction 
2,  for  parallel  geometry,  violates  Woodward-Hoff- 
man  rules.  Hence,  from  the  experimental  evidence, 
we  consider  reaction  I,  summarized  in  equation  6, 
as  the  most  likely  initiation  process.  The  rate  con¬ 
stants  derived  on  this  basis  are  listed  in  Table  2.  Due 
to  space  limitations,  the  conditions  of  these  experi¬ 
ments  are  not  given  but  can  be  requested  from  J.  V. 
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TABLE  1 

Mechanism  Used  for  Fitting  [O]  Profiles  from  H2 


1. 

Ha  +  O2  ^  H  -1-  HOa 

ki  =  fitted 

2. 

Ha  +  Oa  ^  OH  -1-  OH 

ka  =  o’’ 

3. 

H2  -t-  O2  0  4-  H2O 

ks  =  o’’ 

4. 

H2  4-  Kr  -►  H  +  H  4-  Kr 

^4  =  8.86  X 

10““exp( -48321  K/r)° 

5. 

O2  4-  Kr  ^  0  4-  0  4-  Kr 

=  2.66  X 

10“®  exp(  - 59380  K/T)’’ 

6. 

H  4-  O2  ^  OH  4-  0 

he  =  1.62  X 

10““  exp(-7474  K/T)“ 

7. 

HO2  4-  Kr  ^  H  4-  O2  4-  Kr 

kj  =  2.0  X  10“®  exp( -24363  K/T)*' 

8. 

OH  4-  H2  ->  H2O  4-  H 

ks  =  3.56  X 

10  “16  T’  6^  exp(  -  1736  K/Tf 

9. 

0  4-  H2  ^  OH  4-  H 

kg  =  8.44  X 

10““  7^164  exp{  - 3167  K/T)° 

10. 

OH  4-  OH  ^  0  4-  H2O 

kio  =  7.19  X 

10“’ii  T’li  exp(1251  K/T)“ 

11. 

OH  4-  0  O2  4-  H 

kii  =  5.42  X 

10“i3  exp(1112  K/r)° 

12. 

H  4-  H2O  -►  OH  4-  H2 

ki2  =  1.56  X 

10  “16  pi  6^1  exp(  -  9249  K/Tf 

13. 

0  4-  H2O  -►  OH  4-  OH 

ki3  =  7.48  X 

10““  T’l  i  exp{  -  7323  K/r)° 

14. 

H  4-  HO2  ^  H2  4-  O2 

bi4  ~  ki/K^^ 

15. 

H  4-  OH  ->  Ha  4-  0 

ki5  =  3.78  X 

10““  T’161  exp(  -2226  K/r)° 

“All  rate  constants  are  in  cm^  molecule  ^  s 
’’Values  taken  to  be  0;  see  text. 

“Refs.  [23,24], 

•’Ref  [22], 

“Ref.  [25]. 

'"Refs.  [14,26]. 

^Calculated  from  equation  7. 

Michael.  Experimental  corroboration  of  the  exclu¬ 
sive  importance  of  reaction  1  is  not  possible  by  ob¬ 
serving  H  atoms  because  reaction  6  is  effectively  in¬ 
stantaneous  with  the  levels  of  O2  necessary  to 
measure  the  title  reaction.  Experimental  corrobo¬ 
ration  by  observing  OH  is  impractical  since  OH  is  a 
product  in  both  schemes,  also,  the  low  sensitivity  for 
detection  in  comparison  to  ARAS  precludes  design¬ 
ing  a  chemical  isolation  experiment,  that  is,  reactions 
10  and  11  then  become  important. 

Equation  6  can  be  compared  to  earlier  evalua¬ 
tions.  This  equation  gives  values  that  agree  well 
within  a  factor  of  2  with  those  of  Tsang  and  Hamp- 
son  [14]  over  the  experimental  T  range,  1600-2100 
K.  The  values  from  Koike  [16]  range  from  3  to  20 
times  higher,  and  those  from  Skinner  et  al.  [28]  are 
substantially  lower  than  those  calculated  from  equa¬ 
tion  6.  Erom  our  experimental  results,  we  therefore 
conclude  that  reaction  1  is  the  most  probable  pro¬ 
cess  responsible  for  initiation  in  the  H2/O2  system. 
In  the  theoretical  discussion  to  follow,  this  conclu¬ 
sion  will  be  further  assessed  in  terms  of  ab  initio 
determinations  of  potential  energy  surfaces  for  all 
processes.  Subsequently,  conventional  transition 
state  theoretical  (CTST)  estimates  of  the  thermal 
rate  behavior  will  be  presented  and  compared  with 
experiment. 

Theory 

The  electronic  structure  calculations  focus  on  lo¬ 
cating  transition  states  for  reactions  1  to  3.  Hence, 


preliminary  searches  were  made  by  employing  low 
levels  of  theory  which  allow  large  volumes  of  config¬ 
uration  space  to  be  explored,  and  then,  when  a  given 
transition  state  is  located,  the  structure  is  reoptimi¬ 
zed  at  higher  and  more  accurate  levels  of  theory.  All 
calculations  employ  the  Dunning  correlation-con¬ 
sistent  basis  sets  [29-31].  The  preliminary  searches 
were  made  using  the  polarized  valence  double  zeta 
basis  set  (cc-pvdz)  and  the  critical  stationary  points 
were  re-examined  using  the  polarized  quadruple 
zeta  (cc-pvqz)  basis  set.  Searches  were  carried  out 
with  both  coupled  cluster,  CCSD(T)  [32],  and  com¬ 
plete  active  space,  self-consistent  field,  CASSCE 
[33,34],  calculations.  The  CCSD(T)  method  is  the 
most  accurate,  single-reference,  ab  initio  method 
currently  available.  The  CASSCE  calculations  were 
done  to  check  on  the  possibility  of  transition  state 
structures  located  in  regions  not  well  described  by 
single-reference  wave  functions.  In  the  CASSCE 
calculations  reported  here,  all  valence  electrons 
were  taken  as  active  except  the  oxygen  2s  electrons. 
The  active  space  thus  consists  of  8  orbitals  and  10 
electrons  (1512  configuration  state  functions).  All 
calculations  were  carried  out  using  the  MOLPRO 
package  of  codes  [35]. 

A  transition  state  for  reaction  1  was  readily  lo¬ 
cated.  At  the  highest  level  of  theory  used  (i.e., 
CCSD(T)/cc-pvqz),  the  structure  is  as  follows:  Rqo 
=  1.30A,  RoHa  =  I'OSA,  RoHb  =  2.24A,  HaOOZ 
=  107.1°,  and  HbOOZ  =  107.7°.  As  expected  for 
an  endothermic  reaction,  the  transition  state  is  very 
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TABLE  2 

High  Temperature  Rate  Data  for  Hj  +  O2 


T/K 

ki/(cm^  molecule 

2097 

3.0  (-16) 

1927 

9.0  ( - 17) 

1897 

9.0  ( - 17) 

1907 

2.0  (-16) 

1848 

1.0  (-16) 

206.3 

2.5  (-16) 

1981 

1.6  (-16) 

1826 

1.0  (-16) 

1960 

1.0  (-16) 

1819 

3.2  (-17) 

1851 

3.0  ( - 17) 

1891 

1.1  (-16) 

1810 

8.0  ( - 17) 

1679 

5.0  ( - 17) 

1948 

1.5  (-16) 

1974 

1.1  (-16) 

1963 

1.3  (-16) 

2002 

1.4  (-16) 

1881 

8.0  ( - 17) 

1844 

9.0  ( - 17) 

1698 

1.5  (-17) 

1838 

4.2  (-17) 

1819 

2.6  (-17) 

1785 

2.0  ( - 17) 

1662 

1.6  (-17) 

1704 

1.5  (-17) 

1950 

9.0  ( - 17) 

1822 

4.0  ( - 17) 

“late”  with  an  OH  bond  extension  of  only  0.08A, 
relative  to  HO2  and  an  HH  bond  extension  of  0.45A, 
relative  to  H2.  The  calculated,  in-plane,  harmonic 
frequencies  are  1717,  1374,  1169,  321,  and  1825i 
cm“^.  These  can  be  compared  to  calculated  fre¬ 
quencies  for  H2  +  O2  of  4403  and  1600  cm“^  and, 
for  HO2,  of  3672,  1441,  and  1142  cm“^.  For  the 
out-of-plane  bending  mode  of  the  transition  state,  it 
was  not  possible  to  do  calculations  at  the  CCSD(T)/ 
cc-pvqz  level.  For  this  mode  we  use  the  result  of  a 
CCSD(T)/cc-pvdz  calculation,  namely,  681  cm“^. 

The  calculated  barrier,  without  zero  point  correc¬ 
tions,  is  57.8  kcal  mol”^.  For  comparison,  calcula¬ 
tions  with  cc-pvdz  and  cc-pvtz  basis  sets  give  barriers 
of  60.2  and  58.3  kcal  moI“^.  Most  of  the  decrease 
in  barrier  height  with  increasing  basis  set  size  cor¬ 
relates  with  a  decreasing  endothermicity,  that  is,  the 
small  reverse  barrier  is  less  sensitive  to  the  level  of 
theory.  The  calculated  reverse  barriers  are  2.9,  2.6, 
and  2.7  kcal  mol“^  using  cc-pvdz,  cc-pvtz,  and  cc- 
pvqz  basis  sets,  respectively. 

Regarding  the  cc-pvqz  zero  point  corrected  en¬ 
dothermicity  for  reaction  1,  55.4  kcal  mol”^  is  ob¬ 
tained,  and  this  value  is  in  excellent  agreement  with 


recent  measurements  by  Litorja  and  Ruscic  [36], 
55.7  ±  0.8  kcal  mol“^.  This  new  endothermicity 

[36]  is  2.8  kcal  mol“  ^  higher  than  used  by  Janaf  [11], 
and  therefore  the  Janaf  estimated  equilibrium  con¬ 
stants  are  ~1.64  to  100  times  larger  than  implied  by 
Litorja  and  Ruscic,  from  2500  down  to  300  K.  With 
the  harmonic  oscillator  approximation,  for  re¬ 
action  1  can  be  calculated  using  these  experimental 
results  [36],  giving  to  within  ±3% 

Kgq  =  0.2021  7°  exp(- 27656  K/T)  (7) 

between  300  and  2500  K.  can  also  be  calculated 
from  the  present  theory,  and  there  is  only  a  few  per¬ 
cent  difference  between  theory  and  experiment,  cor¬ 
roborating  the  accuracy  of  the  present  ab  initio  re¬ 
sults. 

Extensive  searches  (i.e.,  ~10®  points)  of  both  pla¬ 
nar  and  non-planar  (including  perpendicular)  ge¬ 
ometries  were  also  made  for  the  transition  states  of 
reactions  2  and  3.  The  searches  were  carried  out  at 
both  the  CCSD(T)/cc-pvdz  and  CASSCF/cc-pvdz 
levels.  No  transition  state  for  either  reaction  could 
be  found.  While  the  existence  of  a  transition  state 
connecting  any  two  minima  on  a  multidimensional 
potential  surface  is  guaranteed,  the  existence  of  a 
direct  transition  state  is  not.  Previous  calculations 

[37]  have  shown  the  existence  of  indirect  pathways 
for  both  reactions  2  and  3,  that  is,  H2  +  O2  — >  H 
-I-  HO2  ^  H2O  -I-  O  — >  20H.  In  symmetry-re¬ 
stricted  searches  where  only  geometries  possessing 
a  plane  (or  more)  of  symmetry  are  considered,  high 
lying  saddle  points  can  be  located  for  both  reactions 
2  and  3.  However,  when  the  symmetry  restriction  is 
relaxed  it  is  found  that  these  are  not  true  saddle 
points  because  they  have  two  or  more  imaginary  nor¬ 
mal  mode  frequencies.  The  mode  corresponding  to 
one  of  the  imaginary  frequencies  can  be  identified 
with  reaction  coordinates  for  reactions  2  or  3.  The 
mode  corresponding  to  one  of  the  other  imaginary 
frequency  leads  down  to  the  transition  state  for  1. 
This  suggests  a  potential  energy  landscape  in  which 
the  H2  +  O2  minimum  is  surrounded  by  a  multi¬ 
dimensional  ridge,  the  low  point  of  which  is  the  tran¬ 
sition  state  for  reaction  1.  Thus,  the  picture  that 
emerges  from  these  calculations  is  that  the  lowest 
energy  pathways  for  reactions  2  and  3  are  stepwise 
processes  in  which  the  first  step  is  reaction  1. 

The  final  part  of  this  theoretical  treatment  is  a 
CTST  calculation  of  the  rate  of  reaction  1  using  the 
unadjusted  ab  initio  barrier  height,  structures,  and 
frequencies  given  above.  We  note  that  the  RMS  de¬ 
viation  between  the  calculated  and  observed  fre¬ 
quencies  for  the  reactants  and  products  is  less  than 
30  cm“^,  justifying  the  use  of  unsealed  frequencies. 
These  calculations  were  carried  out  using  the  VAR- 
IFLEX  program  [38],  within  a  rigid-rotor,  harmonic 
oscillator  approximation.  Tunneling  estimates  were 
made  using  both  Wigner  and  Eckart  methods.  The 
Wigner  and  Eckart  tunneling  estimates  are  found  to 
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ELEMENTARY  REACTION  KINETICS 


Eig.  2.  Arrhenius  plot  of  the  data  for  k^.  Open  circles, 
data  from  Table  1.  The  line  is  the  ah  initio  theoretical  cal¬ 
culation  summarized  by  the  three-parameter  expression, 
equation  8.  Closed  square,  data  of  Ref  [13],  and  closed 
circle,  data  of  Refs.  [7-10],  both  transformed  using  equa¬ 
tion  7.  Inset:  open  circles,  data  from  Table  2.  The  line  is 
the  theoretical  calculation  summarized  by  equation  8. 


Eig.  3.  A  simulation  of  OH  profiles  at  three  tempera¬ 
tures  (from  left  to  right,  1927,  1700,  and  1300  K,  respec¬ 
tively),  using  the  mechanism  of  Table  1  with  the  theoretical 
value  for  ki  as  given  by  equation  8.  The  total  density  is 
1.326  X  10^®  with  [Ha]  =  6.524  X  10“  and  [Oa]  =  6.522 
X  10^^  (all  in  molecules  cm“^)  for  all  three  simulations. 
The  simulation  at  1927  K  reproduces  the  experimental  data 
reported  as  figure  4  from  Ref.  [48] . 

be  very  similar,  increasing  the  rate  at  300  K  by  a 
factor  of  3,  while  contributing  only  10%  at  1000  K 
and  less  than  1%  at  2000  K.  Both  the  experimental 
(points)  and  the  theoretical  (line)  results  of  the  pres¬ 
ent  study  are  shown  in  Fig.  2  as  an  Arrhenius  plot. 
The  theoretical  line  in  the  figure  can  be  expressed 
to  within  ±  2%  by  the  three-parameter  equation: 


=  1.228  X  10-1®  y2.4328 

exp  (  —  26926  K/T)  cm®  molecule-®  s-®  (8) 

over  the  temperature  range,  400-2300  K.  Experi¬ 
mental  data  from  previous  work  [7-10,13]  on  the 
reverse  reaction  (  —  1)  have  been  transformed  from 
k_i  to  ki  using  the  revised  K  equation  7.  Note 
that  the  values  transformed  refer  only  to  the  back 
reaction  as  reported  [7-10,13]  and  do  not  include 
the  predominate  process,  H  +  HO2  — >  20H  (esti¬ 
mated  to  be  92  ±  4%  at  room  temperature  [7-10]). 
The  completely  ah  initio  result  summarized  by  equa¬ 
tion  8  is  remarkably  consistent  with  experiment  over 
the  entire  T  range. 

Conclusion 

On  both  experimental  and  theoretical  grounds,  we 
therefore  conclude  that  reactions  2  and  3  can  be 
ignored  and  that  only  reaction  1  initiates  chain 
branching  in  the  H2/O2  system.  The  theoretical  pre¬ 
diction  of  the  rate  behavior  derived  from  first  prin¬ 
ciples  agrees  with  experiment  over  an  extended  T 
range.  We  therefore  recommend  that  the  theoretical 
result,  equation  8,  should  be  exclusively  used  in 
modeling  the  branching  chain  oxidation  of  H2.  The 
remaining  question  to  be  answered  is  then  whether 
this  ki  value  can  supply  sufficient  O  atom  rates  for 
the  branching  chain  oxidation  to  occur.  As  is  well 
known,  ignition-delay-type  experiments  have  tradi¬ 
tionally  monitored  OH  radical  production  [23,39- 
48].  Using  the  mechanism  of  Table  1  along  with 
equation  8,  we  simulated  typical  oxidation  experi¬ 
ments  at  three  temperatures.  These  simulations  are 
shown  in  Fig.  3,  and,  in  all  cases,  chain  branching 
was  clearly  demonstrated.  The  experimental  condi¬ 
tions  for  the  1927  K  simulation  are  from  this  labo¬ 
ratory  [48],  and  it  is  worth  noting  that  the  present 
simulation  reproduces  the  experimental  data  shown 
in  figure  4  from  this  earlier  study.  The  other  two 
simulations  were  carried  out  at  the  same  total  den¬ 
sity,  [H2],  and  [O2]  with  the  only  change  being  in  the 
temperature.  Even  at  1300  K  with  an  incomplete 
mechanism,  rates  of  H  and  HO2  formation  calcu¬ 
lated  from  equation  8  are  large  enough  to  initiate 
branching  chain  oxidation,  although  requiring  a  con¬ 
siderably  longer  induction  time.  Hence,  we  con¬ 
clude  that  reaction  1  is  sufficient  to  initiate  H2  oxi¬ 
dation  by  O2. 
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COMMENTS 

F.  R.  Westmoreland,  University  of  Massachusetts,  USA.  REFERENCE 

What  value  was  used  for  the  heat  of  formation  of  HO2? 

1.  Litorja,  M.,  and  Ruscis,  B.,  ].  Electron  Spectrosc. 
Authors  Reply.  The  value  used  is  from  M.  Litorja  and  97:131-146  (1998). 

B.  Ruscis  [1];  they  report  dHf  Ok  —  (4-0  ±  0.8)  kcal  mole 
(zlHf298K  —  kcal  mole“^). 


